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Abstract 


The  General  Rotorcraft  Aeromechanical  Stability  Pro¬ 
gram  (GRASP)  was  developed  to  analyse  the  steady-state 
and  linearized  dynamic  behavior  of  rotorcraft  in  hovering 
and  axial  flight  conditions.  Because  of  the  nature  of  prob¬ 
lems  GRASP  was  created  to  solve,  the  geometrically  non¬ 
linear  behavior  of  beams  is  one  area  in  which  the  program 
must  perform  well  in  order  to  be  of  any  value.  Numerical 
results  obtained  from  GRASP  are  compared  to  both  static 
and  dynamic  experimental  data  obtained  for  a  cantilever 
beam  undergoing  large  displacements  and  rotations  caused 
by  deformation.  The  correlation  is  excellent  in  all  cases. 

/ 


Nomenclature 


b£  =  basis  vectors  at  the  deformed  beam  tip 
b?*  =  basis  vectors  at  the  beam  root 
e  =  beam  width 

Cij  =  direction  cosine  matrix  relating  the  beam  root 
and  tip 

e  =  error  in  beam  croes-sectional  measurement 
E  as  modulus  of  elasticity 
h  —  geometrical  cross-sectional  property 
/<  =  geometrical  cross-sectional  property 

p  =  projection  of  bf  onto  the  bf  -  bj  plane 
t  =  beam  thickness 
U  =s  strain 'energy 
a  =  material  nonlinearity  coefficient 
0  a=  torsional  deflection  which  was  experimentally 
measured 

k  =  curvature  of  the  beam 
<(>i  =  Rodrigues  parameters 

8  -  load  angle 


Introduction 

The  General  Rotorcraft  Aeromechanical  Stability  Pro¬ 
gram  (GRASP) 1  is  capable  of  treating  the  nonlinear  static 
and  linearized  dynamic  behavior  of  structures  represented 
by  collections  of  rigid-body  and  beam  elements  that  may 
be  connected  in  an  arbitrary  fashion  and  are  permitted  to 
have  large  relative  motions.  GRASP  was  developed  primar¬ 
ily  for  analysis  of  the  steady-state  and  linearized  dynamic 
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behavior  of  rotorcraft  in  the  hovering  and  axial  flow  flight 
conditions. 

GRASP  combines  multibody  and  finite-element  tech¬ 
nology  by  taking  the  strong  points  from  each  area  and  in¬ 
tegrated  them  together  into  a  single,  comprehensive  pack¬ 
age.  GRASP  differs  from  standard  multibody  programs  by 
considering  flexible-body  and  aeroelastic  effects,  including 
simple,  nonlinear,  unsteady  aerodynamics.  GRASP  differs 
from  standard  finite-element  programs  by  allowing  multi¬ 
ple  levels  of  substructures  in  which  the  substructures  can 
move  and/or  rotate  relative  to  others  with  no  small-angle 
approximations.  This  capability  facilitates  the  modeling 
of  rotorcraft  structures,  including  the  rotating/nonrotating 
interface,  and  details  of  the  blade/ root  kinematics  for  var¬ 
ious  rotor  types.  An  overview  of  the  features  of  GRASP 
as  a  program  can  be  found  in  Ref.  1,  while  details  of  the 
analysis  are  given  in  Ref.  2.  The  theoretical  basis  of  the 
analysis  are  addressed  in  companion  papers.3,4 

Because  of  the  nature  of  the  problems  GRASP  was  cre¬ 
ated  to  solve,  the  geometrically  nonlinear  behavior  of  beams 
is  one  area  in  which  the  program  must  perform  well  in  or¬ 
der  to  be  of  value.  The  main  structural  element  in  GRASP 
is  the  aeroelastic  beam,  a  geometrically  nonlinear  beam  el¬ 
ement  based  on  the  kinematics,  internal  and  inertial  forces 
of  Ref.  5  and  the  aerodynamics  of  Ref.  2.  The  element 
kinematics  are  valid  for  small  strains  and  large  rotations, 
but  shear  strains  are  neglected.  The  beam  element  degrees 
of  freedom  include  (as  generalized  coordinates),  rigid-body 
translations  and  rotations  that  are  associated  with  the  stan¬ 
dard  cubic  and  linear  shape  functions.  These  account  for 
the  usual  twelve  degrees  of  freedom  found  in  beam  finite- 
element  analyses.  In  GRASP,  however,  an  arbitrary  num¬ 
ber  of  additional  generalized  coordinates  associated  with 
higher-order  polynomials3  are  also  included,  allowing  the 
element  to  have  a  number  of  degrees  of  freedom  that  is  ar¬ 
bitrarily  more  than  the  usual  twelve,  following  the  so-called 
“p-version"  of  the  finite  element  method.7, 3 

This  paper  is  to  present  numerical  results  from  GRASP  . 
for  comparison  with  static  and  dynamic  experimental  data 
for  large  deflections  of  an  end-loaded  cantilevered  beam. 
(The  experimental  apparatus  and  procedures  are  only 
briefly  described  in  this  paper.  Details  of  the  experiment 
may  be  found  in  Ref.  9  and  10.)  Determination  of  the 
beam  properties  for  input  into  GRASP  is  described  in  de¬ 
tail.  GRASP  results  are  then  presented  along  with  results 
from  previous  analyses  which  are  shown  for  comparison. 
While  the  present  results  do  not  exercise  many  of  the  fea¬ 
tures  and  power  of  GRASP,  they  do  serve  to  validate  much 
of  the  code  dealing  with  the  beam  element's  ability  to  model 
highly  nonlinear 


behavior. 


Experiment 

An  experiment  done  it  Princeton  University8*10  (un¬ 
der  Aeroflightdynamics  Directorate  sponsorship),  was  se¬ 
lected  as  a  test  case  with  which  to  validate  GRASP.  This 
experiment  consisted  of  measuring  the  static  deformation 
and  fundamental  flatwise  and  edgewise  natural  frequencies 
of  a  uniform,  nonrotating,  cantilever  beam  with  a  mass 
attached  to  the  tip  (Fig.  1).  The  beam  was  slender  and 
sufficiently  flexible  to  undergo  large  displacements  (still  at 
small  strains)  due  to  the  presence  of  the  tip  mass.  Beam 
load  angle  and  mass  of  the  tip  weight  were  varied  through¬ 
out  appropriate  ranges. 

The  beam  was  instrumented  with  strain  gages  mounted 
at  the  root  in  the  proper  orientation  to  measure  flatwise 
and  edgewise  natural  frequencies.  The  end  of  the  beam 
was  securely  mounted  in  a  precision  indexing  chuck  that 
provided  a  stable  mount  and  accurate,  repeatable  angular 
settings.  The  static  vertical  and  horizontal  measurements 
were  made  with  a  caliper  scale,  measuring  from  a  flat  table 
with  a  reference  grid  affixed.  The  static  torsional  measure¬ 
ments  were  made  with  the  aid  of  lightweight  reference  rods 
attached  along  the  length  of  the  beam  and  perpendicular 
to  it.  (See  Ref.  10  for  complete  details.) 


though  non-sero  deflections  should  have  been  measureable), 
it  was  necessary  to  assume  that  the  experimental  data  have 
deflections  for  the  no-tip-mass  case  subtracted  out.  This 
latter  point  is  not  explicitly  stated  in  Ref.  9  or  10,  but 
from  the  present  investigation,  appears  to  be  true.  Finally, 
to  allow  for  other  effects  not  present  in  the  simple  elastica 
model,  an  extra  parameter  in  the  form  of  a  material  nonlin¬ 
earity  coefficient  a,  was  introduced.  With  the  assumption 
that  the  beam  is  inextensible,  the  strain  enegy  is  expressed 
as 

U  =  (1) 

where  <c  is  the  curvature  of  the  beam,  /j  and  f4  are  geo¬ 
metrical  cross-sectional  properties  defined  as  follows.  For 
edgewise  deflection: 


r  _  £2 

4  80 

and  for  flatwise  deflection: 


The  beam  was  fabricated  from  7075  aluminum.  The 
mass  density  was  assumed  to  be  2.626  x  10~4lb.-sec.2/in.4. 
The  length  of  the  beam  was  measured  to  be  19.985  in. 
The  thickness  and  width  of  the  beam  were  measured  at 
0.1251  in.  and  0.4999  in.,  respectively.  Assuming  a  gravita¬ 
tional  constant  equal  to  386.089  in./sec.2,  the  mass  per  unit 
length  was  determined  to  be  1.6424  xl0_®  lb.-sec.2/in.2. 
The  mass  momenta  of  inertia  were  2.1420  x  10~*  lb. -sec. 3 
(flatwise)  and  3.4204  xlO-7  lb.-sec.2  (edgewise). 

Determination  of  the  appropriate  values  of  bending 
stiffness  proved  to  be  more  difficult.  Both  static  and  dy¬ 
namic  predicted  behavior  are  sensitive  to  the  value  of  the 
stiffnesses,  therefore  stiffnesses  were  determined  as  accu¬ 
rately  as  possible.  Attempted  inference  of  equivalent  beam 
properties  from  classical  linear  formulas  for  deflection  vs 
load  for  the  two  uncoupled  cases,  load  angles  of  0  and 
90s,  yields  contradictory  information  —  even  when  only 
small  deflections  are  considered.  At  a  load  angle  of  0s  (the 
edgewise-bending  case),  linear  theory  is  too  stiff.  But,  at 
a  load  angle  of  90s  (the  flatwise-bending  case),  linear  the¬ 
ory  is  too  soft.  Similar  contradictory  information  results 
when  attempting  to  correlate  natural  frequencies  of  the  un¬ 
loaded  beam  according  to  linear-beam  theory  with  exper¬ 
iment.  This  suggests  that  there  is  no  one  value  of  E  that 
will  yield  accurate  flatwise  and  edgewise  bending  stiffnesses 
if  the  measured  croee- section  dimensions  are  taken  at  exact 
and  the  theory  is  assumed  to  be  linear.  With  the  failure 
of  linear  theory,  we  turned  to  a  simple,  nonlinear,  planar 
elastica  model.11  Even  then,  a  standard  value  of  the  modu¬ 
lus  of  elasticity  and  measured  cross-sectional  dimensions  in 
a  standard  elastica  model,  only  fair  agreement  is  obtained 
with  planar,  experimental,  static  deflections. 

With  the  aforementioned  problems  in  mind,  a  more  in¬ 
novative  approach  was  called  for.  First,  it  was  determined 
that  only  static  data  should  be  used  because  of  supposed 
accuracy.  Second,  because  transferee  deflections  for  the  un¬ 
coupled,  no-tip-maee  cases  were  recorded  to  be  zero  (even 


'«  =  io 

where  c  and  t  are  the  height  and  thickness  dimensions,  re¬ 
spectively,  of  the  cross-section  of  the  beam.  Considering 
only  the  uncoupled,  static  deflections  (load  angles  of  0  and 
90°),  an  equation  for  the  deflection  of  the  tip  of  the  beam 
was  derived  as  a  function  of  the  beam  bending  stiffness 
and  other  unknowns.  Once  the  equation  for  the  tip  de¬ 
flection  was  derived,  then  a  nonlinear  least  squares  method 
was  used  to  determine  the  best  EIj  and  a  that  fit  the  ex¬ 
perimental  data  for  both  uncoupled  flatwise  and  uncoupled 
edgewise  deflections.  The  value  obtained  for  a  was  ignored 
since  GRASP  does  not  consider  material  nonlinearity.  The 
two  values  of  E  inferred  from  the  bending  stiffnesses  and 
the  cross-section  geometry  were  averaged  and  multiplied  by 
the  cross-sectional  area  to  obtain  the  axial  stiffness.  A  value 
of  Poisson's  ratio  equal  to  0.31  was  assumed  and  the  shear 
modulus,  G,  was  inferred  from  E.  The  following  stiffnesses 
resulted: 

axial  stiffness  =  6.2856  x  10s  lb 
flatwise  stiffness  =  8.4487  x  102  lb-in2 
edgewise  stiffness  =  1.2689  x  104  lb-in2 
torsional  stiffness  =  1.0538  x  10*  lb-in2 

The  ratio  of  the  edgewise  stiffness  to  the  flatwise  stiff¬ 
ness  should  be 

£  =  (4) 

However,  this  ratio,  based  on  the  above  stiffnesses  reported 
above,  turned  out  to  be 


incident  with  the  node  TIP.  Ita  definition  specifies  the  mass 
(5)  and  the  mui  momenta  of  inertiaa  about  all  three  principal 


126.89 

8.4487 


=  15.02 


If  one  aaaumee  that  thia  diacrepancy  ia  due  to  variationa 
in  width  and  thickneaa  along  the  length  of  the  beam,  then 
the  next  question  ia  how  much  variation  would  it  take  to 
cause  this  discrepancy?  Assuming  an  error  e  in  each 
measurement,  and  substituting  in  the  measured  values  for 
e  and  t,  one  obtains  the  ratio  as 


(0.4999  +  «)2 
(0.1251  +  «)2 


=  15.02 


(6) 


yielding  e  =  0.0052  in.,  a  fairly  small  error.  Thus,  the 
inferred  bending  stiffnesses  are  not  unreasonable. 

The  moments  of  inertia  of  the  various  tip  masses  used 
were  estimated  with  some  gross  assumptions  since  these 
values  are  relatively  unimportant.  The  only  properties  of 
the  tip  masses  stated  in  Ref.  9  and  10  were  the  masses. 
There  was  also  a  photograph  depicting  a  tip  mass  with  a 
hollow  cylindrical  shape.  With  this  information  in  mind, 
several  assumptions  were  made:  the  density  of  the  tip  mass 
was  that  of  steel  (0.284  tb/in.s),  the  inner  radius  was  0.375 
in.,  and  finally,  the  length  of  the  tip  mass  was  equal  to  its 
outer  diameter.  With  these  assumptions,  the  moments  of 
inertia  for  tip  masses  of  varying  sise  were  calculated,  (Table 
!)• 


GRASP  Model 


axes. 

Correlation  of  GRASP  Results  With 
Experiment 

GRASP  expresses  static  rotations  in  terms  of  Ro¬ 
drigues  parameters,12’2  so  a  minor  amount  of  postprocess¬ 
ing  is  needed  to  convert  the  GRASP  output  to  the  projected 
angle,  (0),  as  measured  in  Ref.  9,  (Fig.  3).  Consider  the  or¬ 
thogonal  triad  at  the  root  of  the  beam  that  remains  aligned 
with  the  principal  axes  at  the  root.  Introduce  a  dextral 
triad  of  unit  vectors  associated  with  those  axes  denoted  by 
b*  for  i  =  1,2,  and  3.  Now  consider  a  similar  dextral  triad 
at  the  tip  of  the  beam,  denoted  by  bf  for  t  =  1,2,  and  3, 
where  the  deflections  and  rotations  were  measured  in  the 
experiments.  The  relationship  between  the  triads  is  simply 

bf  =  C<,bf  (7) 

where  a  repeated  index  implies  summation.  A  line  along  the 
width  of  the  cross-section  is  then  aligned  with  bf  =  Cub* . 
Now  consider  the  projection  of  bf  in  the  plane  determined 
by  bf  and  bf  denoted  by  p.  The  expression  for  p  can  be 
easily  determined  as 

P  =bf  -  bf  bfb* 

=C„bf+Cl2bf  u 


The  GRASP  model  for  the  Princeton  experiment  is 
depicted  in  Fig.  2.  Subsystem  PRNCTN,  the  model-type 
subsystem  generated  internally  by  GRASP,  represents  the 
complete  structure.  The  first  explicitly  defined  subsystem 
is  CANTBEAM.  The  frame  of  reference  is  defined  to  be  co¬ 
incident  with  the  model  frame  except  for  a  rotation  about 
the  xj  axis  which  is  interpreted  as  the  beam  load  angle.  The 
subsystem  contains  two  structural  nodes  named  ROOT  and 
TIP.  ROOT  is  coincident  with  the  CANTBEAM  subsys¬ 
tem  frame  of  reference,  and  has  all  of  its  degrees  of  freedom 
prescribed  to  xero  (cantilever  beam  boundary  conditions). 
TIP  is  defined  to  be  located  19.985  in.  from  the  frame  along 
the  zj  axis. 

The  first  child  of  CANTBEAM  is  an  aeroelaetie  beam 
element  named  BEAM.  An  aeroelaetie  beam  connectivity 
constraint  associates  the  element's  root  and  tip  nodes  with 
the  nodes  ROOT  and  TIP  in  the  subsystem  CANTBEAM. 
The  definition  of  the  element  includes  specifying  the  or¬ 
ders  of  the  polynomials  used  to  represent  the  displacements. 
The  typical  approach  in  finite  element  programs  would  be 
to  use  several  elements  with  the  transverse  displacements 
approximated  by  cubic  polynomials,  and  the  axial  displace¬ 
ment  and  torsion  approximated  by  linear  polynomials.  In¬ 
stead,  for  this  analysis  we  use  one  element  with  eighth- 
order  polynomials  for  bending  and  sixth-order  polynomials 
for  axial  displacement  and  torsion.  This  yields  a  total  of 
32  element  degrees  of  freedom  (6  of  which  are  constrained 
out  by  the  clamped-end  condition).  Essentially  the  same 
results  are  obtained  when  the  order  of  each  polynomial  is 
reduced  by  one. 

Subsystem  WEIGHT,  the  second  child  of  CANT¬ 
BEAM,  is  a  h  fid- body  mate  element  that  is  defined  to  be  co¬ 


The  angle  measured  in  the  experments  is  the  angle  between 
p  and  bf .  From  Fig.  3,  it  is  clear  that 


0  =  sin-1 


v/l  -  C,3J 

In  terms  of  Rodrigues  parameters* 

lt~  1  +  5±g±g 

c  -*»  +  *** 

l  +  *i**a**i 


(») 


(10) 


where  di<  ds<  and  ds  are  the  Rodrigues  parameters  associ¬ 
ated  with  the  rotation  of  the  tip  node. 

Also,  all  GRASP  deflections  have  the  deflections  for 
no  tip  mass  subtracted  out  before  the  results  are  plotted 
with  the  experimental  data.  All  frequencies  calculated  by 
GRASP  were  converted  from  rad/sec  to  Hs. 

First,  results  are  presented  for  a  1-lb.  tip  mass.  Figure 
4  shows  the  static  deflections  vs  load  angle.  The  GRASP 
correlation  for  flatwise  and  edgewise  is  excellent.  Results 
from  Ref.  13  and  14  are  shown  here  for  comparison,  with 
results  from  Ref.  14  presented  only  for  the  torsional  de¬ 
flections.  Transverse  displacements  from  Ref.  14  were  only 


*  It  should  be  noted  that  the  matrix  of  direction  cosines 
in  this  work  is  the  transpose  of  the  one  in  Ref.  12  and  that 
the  Rodrigues  parameters  used  in  GRASP  differ  from  those 
of  Ref.  12  by  a  factor  of  2. 


available  for  load  angles  of  30  and  40s,  and  therefore,  a 
complete  load  angle  sweep  could  not  be  shown.  It  should 
be  noted,  however,  that  the  transverse  displacement  results 
did  agree  well  both  with  experiment  and  with  GRASP. 

There  were  no  dynamic  results  from  Ref.  14.  The  cal¬ 
culations  presented  in  Ref.  13  are  based  on  the  equations 
of  Ref.  15.  This  analysis  is  restricted  to  moderate  rotations 
caused  by  deformation,  such  that  the  squares  of  the  rota¬ 
tional  components  are  small  compared  to  unity.  Reference 
13  is  not  as  accurate  throughout  the  entire  range  of  load 
angle  as  GRASP.  For  torsional  deflection,  the  GRASP  cal¬ 
culations  cut  right  through  the  middle  of  the  experimental 
scatter.  The  experimental  scatter  here  is  so  large,  however, 
that  it  is  impossible  to  say  which  curve  best  fits  the  data. 
Figure  5  displays  the  flatwise  and  edgewise  frequencies  vs 
load  angle.  The  GRASP  results  are  only  slightly  offset  from 
the  experimental  values,  and  follow  the  trend  exactly.  The 
average  error  is  approximately  0.5%.  Reference  13  does  not 
pick  up  the  trend  for  the  flatwise  frequency,  however  it  does 
follow  the  trend  for  the  edgewise  frequency.  Reference  14 
does  not  consider  the  dynamics. 

The  2-lb.  tip-mass  results  are  presented  next.  Here  the 
torsional  data  have  much  less  scatter  than  in  the  1-lb.  case. 
Again,  GRASP  correlates  excellently  with  the  static  deflec¬ 
tion  as  shown  in  Fig.  6.  Also  again,  Ref.  13  is  close  but 
tends  to  deviate  through  certain  portions  of  the  load  angle 
sweep.  This  deviation  from  the  data  is  large  for  load  an¬ 
gles  above  40s.  Reference  14,  however,  correlates  quite  well 
with  the  static  data.  Figure  7  shows  the  flatwise  and  edge¬ 
wise  frequencies.  The  GRASP  predictions  are  again  slightly 
low  for  both  of  the  frequencies,  but  follow  the  trends  very 
nicety.  Reference  13  while  matching  the  data  fairly  well  at 
0s  load  angle,  strays  from  the  data  at  higher  load  angles. 

Figure  8  presents  the  static  results  from  the  3-lb.  tip- 
mass  case.  All  three  analyses  appear  to  match  the  data 
well  over  the  range  shown.  However,  Ref.  13  results  are 
available  for  only  a  load  angle  up  to  15°,  while  Ref.  14 
calculated  results  up  to  45s.  The  failure  of  Ref.  14  to 
converge  past  15s  is  probably  a  result  of  the  restriction  to 
moderate  rotation  in  that  analysis.  This  is  not  a  problem 
with  GRASP,  since  GRASP  has  no  such  restriction.  In 
Fig.  0  the  flatwise  and  edgewise  frequency  sweeps  are  shown 
for  the  experimental  data  and  the  GRASP  analysis.  The 
correlation  is  excellent. 

For  the  4-lb.  tip-mass  case  there  was  not  a  very  large 
experimental  sweep  available  because  of  the  large  displace¬ 
ments  which  the  beam  underwent.  GRASP  did,  however, 
correlate  very  well  in  the  statics,  Fig.  10,  and  the  dynam¬ 
ics,  Fig.  11.  Reference  14  also  correlates  well  with  the  static 
torsional  deflection  data. 

Figures  12  and  13  summarise  flatwise  and  edgewise 
frequency  vs  tip  mass.  The  lateral  buckling  load  can  be 
inferred  from  Fig.  12.  As  the  tip  mass  increases,  the  flatwise 
frequency  will  tend  toward  sero.  When  lateral  buckling  is 
imminent,  the  frequency  will  drop  to  sero  very  rapidly  as 
indicated  by  both  GRASP  and  Ref.  13.  The  edgewise 
frequency  shown  in  Fig.  13  does  not  reach  sero  since  there 
is  no  buckling  in  this  mode. 


It  should  be  noted  that  the  analysis  of  Ref.  13  suf¬ 
fers  from  being  restricted  to  moderate  rotations.  Reference 
14  does  much  better  than  Ref.  13  in  predicting  the  be¬ 
havior  of  this  configuration  because  equations  used  therein 
are  specialised  for  the  type  of  structure  used  in  the  experi¬ 
ment.  The  equations  of  Ref.  14  are  essentially  identical  to 
those  of  Ref.  13  except  that  certain  terms  of  third  degree 
in  the  unknowns  are  added  to  the  analysis  based  on  the 
observation  that  the  coefficients  of  those  terms  are  large. 
The  size  of  those  coefficients  is  a  function  of  the  ratio  of 
the  stiffnesses  (thus  depending  on  the  cross-section  geom¬ 
etry).  These  added  terms  would  not  be  appropriate  if  the 
cross-section  geometry  were  such  that  the  stiffnesses  were 
of  the  same  order  of  magnitude.  It  is  important  to  note 
that  the  equations  in  the  GRASP  analysis  do  not  require 
that  terms  be  added  or  removed  in  this  manner.  This  is 
an  important  consideration  for  general-purpose  analyses, 
the  equations  for  which  should  not  have  need  of  alteration 
merely  because  of  changes  in  properties. 

Concluding  Remarks 

GRASP  is  a  general-purpose  program  with  both  the 
detail  and  the  generality  to  accurately  model  the  end- 
loaded  cantilever  beam,  as  presented  herein,  as  well  or  bet¬ 
ter  than  the  special-purpose  analyses  in  Ref.  13  and  14. 
Although  this  experiment  demonstrated  significant  nonlin¬ 
ear  behavior  both  statically  and  dynamically,  GRASP  ac¬ 
curately  predicts  the  results.  The  equations  upon  which 
GRASP  is  based  are  not  restricted  as  far  as  the  magnitudes 
of  displacement  or  rotation.  Only  the  strains  are  required 
to  be  small  compared  to  unity.  The  GRASP  analysis  is 
shown  herein  to  be  valid  for  these  types  of  problems.  As 
noted  in  the  introduction,  however,  the  present  validation 
dose  not  exercise  many  of  the  capabilities  of  the  program. 
Also,  as  pointed  out  in  Ref.  1,  the  analysis  does  need  to  be 
extended  to  treat  beams  for  which  shear  deformation  would 
be  important. 

Acknowledgment 

The  second  author  was  supported  by  grant  E-16-697 
from  the  Georgia  Tech  Research  Center. 

Reference* 

'Hodges,  D.  H.,  Hopkins,  A.  S.,  Runs,  D.  L.,  and  Hin- 
nant,  H.  E.,  "Introduction  to  GRASP  -  General  Rotorcraft 
Aeromechanical  Stability  Program  -  A  Modern  Approach 
to  Rotorcraft  Modeling,’'  Proceedings  of  the  42nd  Annual 
Forum  of  the  American  Helicopter  Society,  Washington, 
D.C.,  June,  1986. 

3 Hodges,  D.  H.,  Hinnant,  H.  E.,  Hopkins,  A.  S.,  and  Runs, 
D.  L.,  "General  Rotorcraft  Aeromechanical  Stability  Pro¬ 
gram  (GRASP)  Theoretical  Manual,"  NASA  TM  in  prepa¬ 
ration,  1987. 

’Hopkins,  A.  S.,  and  Likins,  P.  W.,  "Analysis  of  Struc¬ 
tures  with  Rotating,  Flexible  Substructures,"  Proceedings 
of  the  AIAA  Dynamics  Specialists  Conference,  Monterey, 
California,  April,  1987. 


4Hodges,  D.  H.,  Hopkins,  A.  S.,  and  Kuns,  D.  L.,  “Analy¬ 
sis  of  Structural  with  Rotating,  Fiaxibla  Substructures  Ap¬ 
plied  to  Rotorcraft  Aeroelasticity  in  GRASP,”  Proceedings 
of  the  AIAA  Dynamics  Specialists  Conference,  Monterey, 
California,  April,  1987. 

'Hodges,  D.  H.,  “Nonlinear  Equations  for  Dynamics  of 
Pretwisted  Beams  Undergoing  Small  Strains  and  Large  Ro¬ 
tations,”  NASA  TP-2470,  1085. 

'Hodges,  D.  H.,  “Orthogonal  Polynomials  as  Variable- 
Order  Finite  Element  Shape  Functions,”  AIAA  Journal, 
Vol.  21,  No.  5,  May  1983,  pp.  796  -  797. 

TSsabo,  B.  A.,  “Some  Recent  Developments  in  Finite  El¬ 
ement  Analysis,”  Computer *  and  Mathematic *  with  Appli¬ 
cation*,  Vol.  5,  1979,  pp.  99  -  115. 

'Hodges,  D.  H.,  and  Rutkowski,  M.  J.,  “Free- Vibration 
Analysis  of  Rotating  Beams  by  a  Variable-Order  Finite  El¬ 
ement  Method,”  AIAA  Journal,  Vol.  19,  No.  11,  Nov. 
1981,  pp.  1459  -  1466. 

'Dowell,  E.  H.,  and  Traybar,  J.,  “An  Experimental  Study 
of  the  Non-linear  Stiffness  of  a  Rotor  Blade  Undergoing 
Flap,  Lag  and  Twist  Deformations,"  AMS  Report  No.  1194, 
Princeton  University,  Jan.  1975. 


i0Dowell,  E.  H.,  and  Traybar,  J.,  “An  Experimental  Study 
of  the  Non-linear  Stiffness  of  a  Rotor  Blade  Undergoing 
Flap,  Lag  and  Twist  Deformations,”  AMS  Report  No.  1257, 
Princeton  University,  Dec.  1975. 

11  Love,  A.  E.  H.,  A  Treatite  on  the  Mathematical  Theory 
of  Elasticity,  Fourth  Ed.,  Dover  Publications,  New  York, 
1944,  pp.  381  -  426. 

"Kane,  T.  R.,  Likins,  P.  W.,  and  Levinson,  D.  A.,  Space¬ 
craft  Dynamic*,  McGraw-Hill,  New  York,  1983,  pp.  1-90. 

‘'Dowell,  E.  H.,  Traybar,  J.,  and  Hodges,  D.  H.,  “An 
Experimental-Theoretical  Correlation  Study  of  Non-linear 
Bending  and  Torsion  Deformations  of  a  Cantilever  Beam,” 
Journal  of  Sound  and  Vibration,  Vol.  50,  No.  4,  1977, 
pp.  533  -  544. 

"Rosen,  A.,  and  Friedmann,  P.,  “The  Nonlinear  Behavior 
of  Elastic  Slender  Straight  Beams  Undergoing  Small  Strains 
and  Moderate  Rotations,”  Journal  of  Applied  Mechanic*, 
Vol.  46,  March  1979,  pp.  161  -  168. 

"Hodges,  D.  H.,  and  Dowell,  E.  H.,  “Nonlinear  Equations 
of  Motion  for  the  Elastic  Bending  and  Torsion  of  Twisted 
Nonuniform  Rotor  Blades,”  NASA  TN  D-7818,  1974. 


Table  1  Estimated  Inertial  Properties  of  Tip  Mass 


Weight 

-IN 

Lateral  Moments  of  Inertia 
(lb.-in.-eec.') 

Axial  Moment  of  Inertia 
(lb.-in.-sec.*) 

1.0 

1.0822x10-* 

8.2356  x  10“4 

2.0 

3.3676x10”* 

2.6784  xl0“* 

3.0 

6.5673x10"* 

5.3169  xl0“* 

4.0 

1.0561X10-* 

8.6363  xl0“* 

5.0 

1.5276x10“' 

1.2573  x  10“' 

6.0 

2.0658x10“' 

1.7083  x  10“' 

7.0 

2.6670x10“' 

2.2131  x  10“' 

8.0 

3.3278x10“' 

2.7691  x  10“' 

9.0 

4.0457x10“' 

3.3741  x  10“* 

10.0 

4.8185x10“' 

4.0261  x  10“' 

10.42 

5.1589x10“' 

4.3135  x  10“* 

10.46 

5.1919x10“' 

4.3413  x  10“' 

(AEBEAM) 


(RBMASS) 


Fig.  2  Hierarchical  GRASP  model  of  the  Princeton  beam. 
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Fig.  4  GRASP  correlation  with  the  Princeton  experiment 
(1-lb.  tip  hum):  static  deflections. 

a)  flatwise  vs  load  angle. 

b)  edgewise  vs  load  angle. 

c)  torsional  vs  load  angle. 


Fig.  3  Schematic  of  undeformed  and  deformed  beams  show¬ 
ing  bf  (along  the  fixed  axes  at  the  root),  bf  (along  the 
principal  axM  at  the  tip),  p  (the  projection  of  bf  onto 
the  plane  determined  by  bf  and  bf ),  and  0  (the  angle 
meausured  in  the  experiment). 
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Fig.  7  GRASP  correlation  with  the  Princeton  experiment 
(2-lb.  tip  mw):  first  flatwise  and  first  edgewise  frequencies 
vs  load  angle. 

Fig.  8  GRASP  correlation  with  the  Princeton  experiment 
(3-lb.  tip  mass):  static  deflections. 

a)  flatwise  vs  load  angle. 

b)  edgewise  vs  load  angle. 

c)  torsional  vs  load  angle. 
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Fig.  9  GRASP  correlation  with  the  Princeton  experiment 
(3-lb.  tip  mass):  first  flatwise  and  first  edgewise  frequencies 
vs  load  angle. 
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Fi*-  10  GRASP  correlation  with  the  Princeton  experiment 
(4-Ib.  tip  mass):  static  deflections. 

a)  flatwise  vs  load  angle. 

b)  edgewise  vs  load  angle. 

c)  torsional  vs  load  angle. 
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Fig.  12  GRASP  correlation  with  the  Princeton  experiment: 
first  flatwise  frequency  vs  weight  of  tip  ma as. 


Fig.  11  GRASP  correlation  with  the  Princeton  experiment  Fig.  13  GRASP  correlation  with  the  Princeton  experiment: 

(4-lb.  tip  mass):  first  flatwise  and  first  edgewise  frequencies  first  edgewise  frequency  vs  weight  of  tip  mass, 

vs  load  angle. 


